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Abstract 

C**) ' We show that quark orbital angular momentum is directly related to off-forward correlation func- 

tions which include intrinsic transverse momentum corresponding to a derivative with respect to the 
transverse coordinates. Its possible contribution to scattering processes is therefore of higher twist 
and vanishes in the forward limit. The relation of OAM to other twist 2 and 3 distributions known 
in the literature is derived and formalized by an unintegrated sum rule. 
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1 Introduction 

It is well known that orbital angular momentum of quarks and gluons can give an important contribution 
to the total spin of the proton, according to the sumrule 

i = ~Aq + L q + Ag + L g , (1) 

where Ag, L q , Ag and L g are the first moments of the corresponding quark spin Aq(x), quark OAM L q (x), 
gluon spin Ag(x) and gluon OAM distribution L g (x). Further on one knows since long how to measure 
polarized distribution functions and much work was invested into reliable estimates of their leading 
moments. Still, despite enormous theoretical and experimental efforts only Ag is known to be rather 
small, and the largest value given in the literature is of the order of Aq(Q 2 ~ lGeV 2 ) ~ 0.30 in 

the MS scheme, while our knowledge of Ag is extremely limited 0] 03 ■ Finally, for a long time it was not 
known at all how to access the quark and gluon OAM pieces experimentally With the advent of studies 
of off-forward scattering processes like DVCS, L q and L g came for the first time within reach, although 
not directly, but via the sumrule |Sj 

~ J dxx(H q>g (x,0 + E q , g (x,0) = J q , g , (2) 
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where H(x,£) and E(x,^) are off- forward distributions in the limit of vanishing momentum transfer 
squared A 2 — > 0. This sumrule allows to extract e.g. the quark OAM contribution, assuming knowledge 
of the quark spin Aq, via 

L q = J q - \A q . (3) 

While an experimental determination of H qg (x, £, A 2 ) and E q:9 (x, £, A 2 ) with an accuracy which allows 
to estimate J q and J g with interesting precision is still a long way to go, eq. was used recently to 
calculate J q on the lattice 

In ref.||] a sum rule has been presented relating quark OAM to the second moment of the twist 3 
off- forward distribution Gs(x) (now G^Oc))- Our notation is that GPDs depending only on x denote 
the GPDs taken in the forward limit, e.g. Gi(x) — Gi(x,£ = 0, A 2 = 0). Including the new Gi(x) 
distribution this sumrule reads 

J dxx [G 2 {x) - 2G 4 (x)] = -L q . (4) 

It is useful to notice that the identification eq.(@J has been made using the sumrule eq.©- On the level 
of distributions, Hoodbhoy et al. showed that the sumrule eq. J5J) is valid for higher moments as well 
and therefore in the forward limit 

±x(q(x)+E(x))=J q (x), (5) 

where it has been used that H{x) is equal to the unpolarized quark distribution q{x). Furthermore 
Hoodbhoy et al. defined quark spin and OAM distributions (in terms of inverse Melin-transformed 
higher moments) which are evidently interrelated by 

L q (x) = J q (x) - ^Aq(x). (6) 

Combining cq.JSJ) and eq.©, it is possible to determine the quark OAM distribution directly from mea- 
surable quantities using the sumrule 

L q (x) = l -x (q(x) + E{x)) - l -Aq(x). (7) 

In this study we show that the parton model definition of the quark OAM distribution can be related to a 
certain off-forward matrix element as soon as intrinsic transverse momenta are taken into account. Using 
the equation of motion, we derive an unintegrated sum rule in the Wandzura-Wilczek approximation 
which relates L q (x) to some of the known twist 2 and 3 off- forward distributions. 



2 Quark OAM, intrinsic transverse momentum and off-forward 
correlators 

2.1 Rewriting Quark OAM 

We begin by recalling the definition of the quark OAM distribution in the light cone gauge according to 

m 

h « {x) = J dx e 2 Mj-d^x) ' (8) 

with = \/2~f~^i + %p and where we have replaced the residual gauge covariant derivative by the partial 
derivative, T> = d — igA — > d. If the boundary conditions can be fixed so that the A^-fields vanish at 
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infinity, then the residual gauge field A would be exactly equal to zero. However in the light cone gauge, 
non- vanishing gluon fields at the boundary can give rise to topological effects, as is discussed in How 
this could effect our calculation lies beyond the scope of this presentation. First we observe that one 
easily runs into problems due to the explicit factors and the integration of x± in case that one takes the 
definition eq.JBJ) literally, see below. Our plan is now to reformulate eq.(JHJ) in order to circumvent this 
from the beginning by incorporating an additional transverse vector which gives us some handle on the 
transverse direction. For this purpose we introduce the function /l 5 (x, Aj_) defined as 

f Lq (x,A ± ) = J dX ~^ X± (P'\ i/>l(x ± )i (x 1 2 - x 2 d,) il>+(x ± + x~) \P) , (9) 

which differs from eq.JSJ) mainly in that the matrix element is slightly off- forward in the transverse 
direction, i.e. P' = P + Aj_/2, P = P — Aj_/2. The function Jl (x, Aj_) is well defined when acting on 
a test function T(AjJ, and we therefore define the quark OAM as 

/L « (a ° " (2^(0) / ^ 2 ^T(A x )/ ig (x,A x ). (10) 

We see that the original definition eq.© can be reproduced by choosing T(Aj_) = S 2 (A±). Shifting the 
fields and writing the factor xj_ as derivative with respect to Aj_, we find after partial integration 

ejkOA ±J {T(A ± ){P'\^ + (0 ± )d, ±k Mx ± + x-)l^ \P)}. (11) 

where the antisymmetric ejk — 1 for j = 1, k = 2 and for j = k. Since we get no contribution in case 
that the derivative acts on the testfunction (e.g. if T(Aj_) is symmetric in Aj_), we end up with 

KW = j ^ x ^e lk d A±1 {(P'\^l(0^d x± Mx^+x-)\ x±=o \P)} A ^ o . (12) 

The choice T(Aj_) = 5 2 (A±) would result in an integral over the square of the delta function, which is 
not well defined. This is an indication for the difficulties associated with a naive use of (jHJ. 

Another way to circumvent these potential problems is to introduce wave packets \<f>), e.g. similar to 
the discussion in |14|. We then would start with the forward distribution /(x) and make the integration 
over the transverse coordinate x± as in eq. JSJ more explicit by writing 

f{x) = J d 2 x±g{x, x x ) = J d 2 Xl _ (P\ x x 6(x, x x ) \P) ■ (13) 

The function g(x,x±) can now be treated in exactly the same way as the impact parameter dependent 
distribution in |14| . and by replacing the momentum eigenstates in ea. (ll3H by wave packets, we find 

/(*) = 7^1 d 2 A^6 2 (A ± ) J d 2 P ± cb*(P^ + ^MP± ^) 

id A± (P± + ^-\d(x,x ± = 0)\P ± -^-) (14) 

where N is normalization of the wave packet. This has to be compared with the corresponding result 
using the testfunction (see the structure of eq. JTJ ), 

/(*) = i / d 2 A ±< 5 2 (A ± )T(A ± )za Ai (P ± + ^\6(x,x ± = 0)\P ± - (15) 



3 



From cqs. J2| and (JTSJ we find that T(A ± ) plays a role similar to / d 2 P^(j)*{P 1 _ + %-)<KPi_ - 4^). 

Orbital angular momentum distribution functions similar to eq. have been used in refs.|15| and 
[Tt)| to calculate their respective evolution equations. A study of the evolution of off-forward correlators 
and the sumrule eq.QJ for dressed quark states in light-front pQCD can be found in |T7) . 



2.2 Quark OAM in terms of proton wave functions 

We now reexamine the above considerations using proton wave functions. Starting from the decomposition 
of the proton state \P) in terms of proton wave functions JS| (f° r the overlap representation of GPDs 
see also (15) 1 



n 



cPki i dxi 



{xiP + , XiP + + fcj_L, Aj} i=1 n ) fcix, A;, A), 



(16) 



where \I/ n (a;i, fcij_, A.;, A) is the n-particle fock state wave function of the proton. The normalization is 
given by 

(p + ,p±,X'\ k+,k ± ,X) = 2 P + (2Tr) 3 S xx ,5( P + - k+)5W( P± - k ± ). 



We choose 



p + = p o + pa^ p- = pQ _ p 3 ; 

P' = P + A/2,P = P- A/2,(1 + £)P'+ = (l-0^ + , 



(17) 



but work most of the time in the limit £ = 0. 
Inserting ea. l|16f) in eq.|0 we end up with 



/£ f (*.ax) = ^EE/ 



n dx l d 2 q l _ 
„i 2 ( 2 ^) 3 



"(f>)/' 



£a 9a), 



(18) 



where 



q a ± = q a ± + (1 - £a)Aj_, g^o-L = 9a - XiAj.. 

We see immediately that the rhs of ea. (|18fl vanishes if we take the naive limit Aj_ —> 0, which would 
reproduce the definition eq.JSJ). Using instead the definition ea. (|10fl . rewriting x± as a derivative and 
performing a partial integration gives 



i i 



T(0) 2 



^EE 



n 



dxid 2 qix 



5 (x a - x) 8 I 1 - ; 



1 2(2^)3 

5 {2) (E^l ({?» x ^Ala^^i^a)^^,^)) | Aj 

+ (xi, aa)*„(asi, gii) {q a x d A } 3 T(A±)) \ A±=0 } , 



(19) 
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where we have dropped all terms oc Aj_. Since the second term in (|19|l vanishes, we end up with 



fZ F W 



Jee 



dx l d 2 q i ± 
i\ 2(2^)3 



n 



6 (x a - x) S I 1 - ^2 ■ 



{x l ,q l ±_)^ n {x ll qi±)) |a_l=0 



(20) 



which is the quark OAM in terms of proton wave functions. This can be written as 

fL q F (*) = J ^ x ^e ±jk d A±j [{P'\^ + (x x )d x ^ + (x ± + x-)\P)) | Ax ,^=o 
which is just equal to eci. ((T2)) . 

2.3 Off- forward correlators and intrinsic transverse momenta 

First we introduce the slightly more general kinematics 

n 2 = p 2 = 0, n ■ p = 1, 

M 2 1 

p = (i -Op + (i + 0-2~n+2 A ^> 

M 2 1 

P = (l + 0P+(l-0-2-»"2 A -L. 



M 2 = M 2 
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(21) 



:A 2 



(In terms of light cone coordinates one can choose 



+ 



P 

p~ = n + — 0.) 
Now let us consider the generic off-forward correlator 



2 
P+ 



dX^± £lx ^ lk± . x± (p/| _^ )w( A nj |, |p) 



(22) 



In order to get a gauge invariant correlator, we included a link operator hi which runs in particular along 
the transverse direction connecting the points (— A/2 n, —x±/2) and (A/2 n, x±/2). For recent discussions 
of these transverse gauge links and their implications see e.g. |2 11122) . Neglecting the intrinsic transverse 
momentum k± in a given hard scattering amplitude involving ca. (|22|l allows for a direct integration of 
eq. (|22|l over k± , and in this case we end up with an expression which is local in the transverse direction, 
x± = 0, 



rfA 
2tt 



A 



A 



(23) 



Such correlations functions have been parametrized in terms of twist 2 and 3 distributions, see e.g.[S||§]. 
Taking into account terms linear in the intrinsic transverse momentum k 1,± and performing the integration 
leads to a correlator with one derivative [23j 



x±=0 ' 



(24) 
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where d = 1/2(9 — d). The partial derivative, when acting on the link operator U in ea. 1)22(1 leads 
to terms with explicit transverse gluon operators A±. These additional contributions will be neglected 
in our approximation. In any case, there are no transverse links left in eqs. since the relevant 

operators are local in transverse direction. For the following it is important to observe that the correlator 
eq. (l24|l must be proportional to the components of the only remaining transverse vector, the momentum 
transfer Aj_. We concentrate now on the part which is proportional to the combination e" ±rT± A a± . Then 
it is possible to parametrize eq.J23J by 

J _ e ^(P'|^(--n,-^) jhd»^(-n,^)\P) x±=0 = -^—U(P\S') jh l5 U(P, S)L q (x, £, A 2 ) 

+ (ocA 1 ^ ) + («£) (25) 

where we have already multiplied with n M . The last two terms on the rhs of ea. (|25|l indicate the presence 
of contributions proportional to A v± and £ which vanish after taking the derivative e Vx<7± dV- (see below) 
and the forward limit. Of course one can write down other terms for the parametrization in ea. 1251) which 
only implicitely give rise to a factor e" ±a± A a± . Using (generalized) Gordon identities (see discussion 
below eq.(j2Bl)) these contributions can be, however, reduced to the term on the rhs of ea. (|25(l . up to the 
indicated terms which are irrelevant for our calculation. We now rewrite the lhs and obtain 

/ ^e <£ ^(-+0 Qa"- + (P'| Vl(0)V+(^ + x-) \P) | Xi=0 (26) 

The function L q in ea. ()25fl can then be extracted by taking the derivative £„ iffl 3J" of both sides. For 
the rhs of 1125H we get in the forward limit 2L q {x) = 2L q (x, £ = 0, A 2 = 0), and the lhs is given by 

/ ip^^Wx^ ((P'l ^(0)^+(*l + x-) \ x±=0 (27) 

Comparing eq.J2U and ea. (|27J) . we see that in the forward limit we have the identification 

L q (x) = f Lq (x). (28) 
Thus the off-forward correlator l|24|) is directly related to the quark OAM distribution of the proton. 

2.4 Relation to other twist 2 and 3 off-forward distributions 

Since the correlator ea. l)24|l involves a transverse derivative, it corresponds to (kinematical) twist 3. 
Similar to the investigations in [201 123| we apply now the identity 

^2)7^^1) = -^"^2)7^75^1), (29) 

£0123 = +1, (30) 

coming from the equations of motion, where d = 1/2(9 + d) and where [•••] stays for the anti- 
symmetric combination of the indices, to the correlator ea. 1(25(1 . We contract again with n M and choose 
v = v±. This leads to the following relation 

/ ^ (P'\^n,-^) 0Z^ n ,^)\P) x±=o 

= / ^e^{P>\j(-±n,-^)^n. ^^ n ,^)\P) x ^ 

_ ie n^ J ^ {J y| ^^^^ ^ | p) ^ =Q . (31) 
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Using 



the third line in eq. Ij31|l can be rewritten 

dX 



X x± 



°J a± n ■ d 



X x± . 



(32) 



2tt 
dX 



e «A (p /| ^ - } jha^^n, -±) \P) 



ii=0 



+« wiaiP y ^e" A (P'| ft-^n, -^) 7 a x n ■ £ 7 5^(^, ^) |P) Xi=0 . (33) 

Using partial integration (when possible) and translations, one can show that the three different deriva- 
tives occurring in eq. l|3T|l and eq. can be substituted by 



n ■ d 


■* —ix, 






n ■ d 


i 


-> —n ■ 
2 


. m 





(34) 



2 
dX 
27 



X 



X x± , 



(P'| ^(--n, —±) M^-n, |P) Xi= 



dA 



2 ' 2 
A x± 



2 ' 2 
A x 



(35) 



The correlators in eq. 1(35(1 can be completely parametrized in terms of the known twist 2 and 3 off- forward 
distributions in the WW-approximation, see ^UJ- This gives 

U(P', S') {-ix [(H + E + G 2 ) 7 ^ - G 4 ie^ np A a± /i l5 ] 
. i 



~ n • A 



(H + G 2 ) la±l5 + (E + G 1 )^ l5 



+G 3 A a± A7s] } U(P, S) + (oc A VJ -) 



(36) 



where all GPDs are functions of x,£ and A, and where we do not show the terms proportional to A v ^ , 
because they are not directly related to the L g -term in ea. <(25|) . Using some Gordon-identities [21] we 
have with our kinematics and in the limit £ — > 

e vx0xnp A 0± U(P',S') jh-f 5 U{P,S) = -2iU(P',S')^U{P,S). 

Taking this together with ea. 1)25(1 we end up with the following sumrule between the distribution L q (x) 
from and the distribution functions from eq. ((36(1 in the forward limit, 



L q {x) = x (q{x) + E(x) + G 2 (x) - 2Gi(x)) - Aq(x). 



(37) 



where L q (x) is, according to the upper analysis, identified as (forward) quark OAM in the parton model. 
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2.5 Comparing with the integrated G2-sumrule 

For reasons of comparison let us now recalculate the integrated sumrule. Following ref. |25| (eq.50), we 
have due to equations of motions 

/ dxx J S e " A ( p/ l^-^-T )r ^ ( 5 n 'T )|p) -=° 

= U(P' , S')^ U± U(P, S)^- I dx [x(H(x) + E(x)) + Aq(x)] + cc A v± . (38) 



2 . 

According to the parametrization in 10], this is equal to 

U(P',S')f Vx U(P,S) J dxx [H(x)+E(x)+G 2 (x)} 
-ie v± P ±np A 0x U(P',S') jh l5 U(P,S) J dxxGi(x) {+ cc A v± ) , 
Taking both equations together we get 

U(P',S')Y X U(P,S) J dx -~x(H(x)+E(x))-xG 2 (x) + ^Aq(x) 
= -ie" ±(3±np At3 ± U(P',S') jhj 5 U(P : S) J dxxGi(x). (39) 
which leads to the sumrule 



dx 



-x{H{x) + E(x)) - xG 2 {x) + 2xG 4 {x) + ^Aq(x) 



0, 



respectively 

dxx [G 2 (x) - 2G 4 {x)} = —J q + ^Ag = -L q . (40) 
On the other hand, 

J dxxG 4 (x) = 0, 

giving the known E3 result 



J dxxG 2 {x) = —J q + -Aq = -L q . 



This sumrule is gauge-invariant, and the possible contributions including explicit gluon operators drop 
out for the second moment over x. Furthermore, integrating 1)37(1 over x and using Ji's sumrule J2J) 6 
we find 



L q = j dxx (q(x) + E(x) + G 2 (x) - 2G 4 (x)) - j dxAq(x) 
= 2J q -Aq + J dxxG 2 {x) 

= 2L q + J dxxG 2 (x), (41) 

which is perfectly consistent with the sumrule eq. (|40|l and confirms our identification eq. (|28|l on the level 
of integrated distributions. 
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3 Conclusions 



The quark OAM distributions in eq. Q and eq. i|37|) as well as the definitions in |15l ITC) all coincide 
when contributions are neglected which contain explicit transverse gluon operators A±. This has been 
discussed in JT], see eqs.(20)-(24) and the paragraph below eq.(28) therein. Dropping the v4^-terms, we 
can therefore combine eq. JjJ and ea. (|37fl to get the interesting and simple relation 

L q {x) = -x [G 2 (x) - 2G 4 (x)] , (42) 

which is obviously a generalization of the integrated sumrule eq. ()40(l . 

In summary we have shown that in the framework of the WW-approximation the quark orbital angular 
momentum distribution is directly related to the twist-3 GPDs G2(x) and Gn(x), taken in the forward 
limit, in form of the sumrule ea. (|42|) . Our results represent only a small step towards solving the notorious 
problem of a direct measurement of the quark OAM contribution to the nucleon spin. At least it leads 
to a new and nice interpretation of the above mentioned twist-3 GPDs. It has to be seen if the sumrule 
ea. (|42[) or a similar expression holds outside the WW-approximation. The main obstacle in this regard 
will be the use of an accurately defined gauge invariant OAM distribution. 
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